TAUP-2554-95 



Generalized Boltzmann Equation 
in a Manifestly Covariant 
: Relativistic Statistical Mechanics 

■ L. Burakovsky* and L.P. Horwitz' 

School of Physics and Astronomy 

CN) ■ Raymond and Beverly Sackler Faculty of Exact Sciences 

<N : Tel- Aviv University, Tel- Aviv 69978, Israel 

oo 
O 
in 
P> 

Abstract 

We consider the relativistic statistical mechanics of an ensemble of N events 
with motion in space-time parametrized by an invariant "historical time" r. 

■ We generalize the approach of Yang and Yao, based on the Wigner distribution 
functions and the Bogoliubov hypotheses, to find the approximate dynamical 
equation for the kinetic state of any nonequilibrium system to the relativistic 
case, and obtain a manifestly covariant Boltzmann-type equation which is a 



relativistic generalization of the Boltzmann-Uehling-Uhlenbeck (BUU) equa- 
tion for indistinguishable particles. This equation is then used to prove the 
-H-theorem for evolution in r. In the equilibrium limit, the covariant forms of 
the standard statistical mechanical distributions are obtained. We introduce 
two-body interactions by means of the direct action potential V(q), where q is 
an invariant distance in the Minkowski space-time. The two-body correlations 
are taken to have the support in a relative 0(2, l)-invariant subregion of the 
full spacelike region. The expressions for the energy density and pressure are 
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obtained and shown to have the same forms (in terms of an invariant distance 
parameter) as those of the nonrelativistic theory and to provide the correct 
nonrelativistic limit. 

Key words: special relativity, relativistic Boltzmann equation, relativistic Maxwell- 
Boltzmann/Bose-Einstein/Fermi-Dirac, mass distribution 
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1 Introduction 

This paper continues a series of works on relativistic kinetic theory of an iV-body 
system within the framework of a manifestly covariant mechanics |§, both 

for the classical theory and the corresponding relativistic quantum theory. In this 
framework, for the classical case, the covariant dynamical evolution of a system 
of N particles is governed by equations of motion that are of the form of Hamil- 
ton equations for the motion of N events which generate the particle space-time 
trajectories (world lines). These events are considered as the fundamental dynam- 
ical objects of the theory and characterized by their positions g M = (ct, q) and 
energy-momenta p^ = (E/c,p) in an 8iV-dimensional phase space. The motion is 
parametrized by a continuous Poincare-invariant parameter r || called the "historical 
time" . For the quantum case, the covariant dynamical evolution of N particles is gov- 
erned by a generalized Schrodinger equation for the wave function Vv(<?i, <?2, Qn) £ 
L 2 (R 4N ), with measure dq\dq2 ■ ■ ■ dq^ = d 4N q, describing the distribution of events 
qi = gf , /i = 0, 1, 2, 3; i = 1, 2, . . . , N. The collection of events (called "concatena- 
tion" ||) along each world line corresponds to a particle in the usual sense; e.g., 
the Maxwell conserved current is an integral over the history of the charged event 
10|| . Hence the evolution of the state of the iV-event system describes a posteriori 



the history in space and time of an iV-particle system. 

The evolution of the system is assumed to be governed by Hamiltonian-type equa- 
tions with a Lorentz-invariant scalar function, the relativistic dynamical function of 
the variables (qi,Pi) specifying the state of each particle i. In the simplest case of a 
free particle, for which the world line is generated by a free event, the relativistic 
dynamical function (generalized Hamiltonian) is 

2M ' 

where we use the metric g^ v = (— , +, +, +), and M is a given fixed parameter (an 
intrinsic property of the event), with the dimension of mass. 
The Hamilton equations 

dq 11 _ dK dp" dK 
dr dp^ dr dq^ 
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M 



yield, in this case, 

~d7 ~ m' ~dT ~ 

Eliminating dr, one finds 

<H £ 2 

dt £° ' 

as required for the motion of a free relativistic particle. It then follows that, for a 
free motion, the proper time interval squared, divided by dr 2 , is 

dxfdq^ _ p^Pf, 
dr dr ~ M 2 

For 

corresponding to the "mass-shell" value 

P % = -mV, 

it follows that 

c 2 dt 2 - dq 2 = c 2 dr 2 . 
In the more general case in which 

K = K + V, 

where V is, for example, a function of q, p 2 = p^p^ may vary from point to point 
along the trajectory. Hence, in general, the proper time interval does not correspond 
to dr. 

For a system of N interacting events (and hence, particles) one takes 

K = Y:^ + V( qi ,q 2 ,...,q N ), (1.1) 

where all of the events are put, for simplicity, to have equal mass parameters, and we 
write qi, for brevity, for the four- vector. The Hamilton equations are 

d$ _dK _pt 
dr dp ifl W 

M = _^ = _^l. (1.2) 
dr dq itJ , dq^' 

These equations are precisely of the same form as those of nonrelativistic Hamilton 
point mechanics, but in a space of 8N dimensions instead of 6iV. The fundamental 
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theorems of mechanics, such as the Liouville theorem [0 , the theory of canonical trans- 
formations and Hamilton- Jacobi theory, follow in the same way, with the manifold of 
space-time replacing that of space, and energy-momentum replacing the momentum. 
It is fundamental to this structure that there is a single universal evolution parameter 
r which plays the role of the Galilean time. 

In the quantum theory, the generalized Schrodinger equation 

d 

^0-Vv(?l, 92, • • • , Qn) = Kil) T {q u q 2 , . . . , q N ), (1.3) 

with, for example, a K of the form (1.1), describes the evolution of the iV-body wave 
function Vv(<?i, <?2, • • • , Qn)- To illustrate the meaning of this wave function, consider 
the case of a single free event. In this case, (1.3) has the formal solution 

VV(?) = (e- tKoT/h ^o)(q) 

for the evolution of the free wave packet. Let us represent ip T (q) by its Fourier 
transform, in the energy-momentum space: 

Ml) = T^U / d 4 pe-^ 2Mn e^ h Mp), 

where p 2 = p^p^, p ■ q = p^q^, and ipo(p) corresponds to the initial state. Applying 
the Ehrenfest arguments of stationary phase to obtain the principal contribution to 
ip T (q) for a wave packet centered at p^ we find 

Qc M 

consistent with the classical equations (1.2). Therefore, the central peak of the wave 
packet moves along the classical trajectory of an event, i.e., the classical world line. 

The wave functions have a local interpretation, i.e., \^ T {.q)\ 2 d A q is the probability 
to find an event at the space-time point q^ in space-time volume d 4 q. Localization in 
space, as well as in time, can be shown by applying arguments given in ref. JTT| . 



Horwitz, Schieve and Piron |J have constructed equilibrium classical and quantum 
Gibbs ensembles. They found that the grand partition function in the rest frame of 
the system is given by 



\nZV},V,n,» K ) = J ^e-^w, p = __. (1.4) 

In addition to the usual chemical potential \x in the grand canonical ensemble, there 
is a new potential fix corresponding to the mass degree of freedom of relativistic 
systems. 

Horwitz, Shashoua and Schieve || have shown that in the framework of the man- 
ifestly covariant mechanics which we discuss here, covariant Weyl transforms exist 



,2 



4 



for observables, and therefore covariant relativistic Wigner functions [|T^] can be con- 
structed. In this way they derived a manifestly covariant relativistic generalization 
of the BBGKY hierarchy for the s-particle relativistic Wigner functions. By approxi- 
mating the effect of correlation of second and higher order by two-body collision terms 
(using the cross-sections defined in ref. as in the usual nonrelativistic Boltzmann 

theory, they obtained a manifestly covariant Boltzmann equation (for non-identical 
events). This equation was used to prove the if -theorem for evolution in r. In the 
equilibrium limit, a covariant form of the Maxwell-Boltzmann distribution, 

f(°\q,p) = e A (i)(p-P°)\ (1.5) 

was obtained. Since this distribution is the distribution of the 4-momenta of the 
events, m 2 = —p 2 = —p^p^ is a random variable in a relativistic ensemble. In order 
to obtain a simple analytic result the authors restricted themselves to a narrow mass 
shell p 2 = —m 2 = —M 2 . The results obtained in this approximation are in agreement 



with the well-known results of Synge for an on-shell relativistic kinetic theory. 

In ref. @ the equilibrium Maxwell-Boltzmann distribution (1.5) was considered 
for the whole range of m, to obtain the corresponding equilibrium relativistic mass 
distribution. Its low-temperature and nonrelativistic limits were investigated and 
shown to yield results in agreement with nonrelativistic statistical mechanics 0[. 

In the present paper we study the case of indistinguishable events. In contrast to 
the approach of Horwitz, Shashoua and Schieve ||, we choose another approach ini- 



tiated by Yang and Yao |L5| in the nonrelativistic case, which is based on the Wigner 



distribution functions and the Bogoliubov hypotheses to find approximate dynamical 



equation for the kinetic state of any nonequilibrium system [16]. Kinetic equation that 
we obtain, which represents a relativistic generalization of the Boltzmann-Uehling- 
Uhlenbeck (BUU) equation [I7| for indistinguishable particles, and can be easily gen- 



eralized to include the non-identical well. The generalized Boltzmann equation 

obtained in this way is then used to prove the if -theorem for evolution in r. In the 
equilibrium limit, the covariant forms of the Bose-Einstein/Fermi-Dirac/Maxwell- 
Boltzmann distributions are obtained, which, as considered for the whole range of m, 
provide the corresponding equilibrium relativistic mass distributions. The relativistic 
mass distributions are studied in the identical particle case in [|J, and their possible 
consequences for high energy physics and cosmology are considered, respectively, in 



and 1 



We introduce two-body interactions by taking the support of mutual correlations 
for any two events to be in a relative 0(2, l)-invariant subregion of the full spacelike 
region, as done in the solution of the two-body bound state problem fTS, ED[ , and for 



the extraction of the partial wave expansion from the relativistic scattering amplitude 



2jJ . We then calculate the expressions for the energy density and pressure of an 
interacting gas, and show that they have the same form (in terms of an invariant 
distance parameter) as those of the nonrelativistic theory and provide the correct 
nonrelativistic limit. 
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2 Relativistic TV-body system 

The evolution in r of an iV-body system is determined by the Liouville-von Neumann 
equation for the iV-body density matrix p (we use the system of units in which 
H = c = ks = 1, unless other units are specified): 

i% = [K,p\, (2.1) 
where K is the total iV-body Hamiltonian, here taken to be 

N N 

K = Y,K? ] + £ Vij, (2.2) 

i=l l=i<j 

where 

i 2M 

and 



Vi,j = V(qt- qj) , qi-qj = \J (q? - qf) (q^ - ty>) 

is a two-body interaction potential. In order to obtain the BBGKY hierarchy, one 
introduces the (n)-body density matrices, as follows: 

N\ 

Tr (l>2,..,n)PlX..,n = ( N _ n)V ( 2 ' 4 ) 



Pll...,n = ( N ^ n]] Tr (n+l,.,N)P, (2.3) 



and, by taking the appropriate traces in Eq. (2.1), obtains p2 



n n 

= E^°Ut.,J + E^UlJ 

^ r j=l l=i<j 



+ Tr {n+1) J2[Vi,n + i,P?C,n + i}- (2-5) 

i=l 

This set of equations is equivalent to (2.1). 

In what follows, we shall use the simplified notation: = pf\ p iy j = p?j , etc., so 
that the latter equation can be rewritten as 

q n n n 

^ = E [#* . Pn] + E Wij> Pn] + Tr (n+1) £ Pn+l]- (2.6) 

OT i=l l=t<j i=l 

It is convenient to introduce directly the symmetry requirements on the function p n 
by means of 

pn = S n F n , (2.7) 
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where S n is a symmetrization/aritisymmetrization operator defined by 



i-1 



S n = H l±£Py > . (2.8) 

i=2 V j=l / 

Here Py denotes the permutation operator. Since S n satisfies the relation 

S n +i = S n ^1 ± ^ fj.n+l^ (2-9) 

and commutes with the operators Ki and Vjj, one can substitute (2.7) into (2.6) and 
obtain the equation 

ap n n n 

ijT = Et^i^n] + E [^,i^n] + Tr (n+1) ^[y i>n+1 ,F n+1 ] 

^ r i=l 1=«<J i=l 

n n 

± T ^(n+1)E [^n+l,E P ^+l^+l]- (2-10) 
i=l i=l 

Now we introduce the Wigner distribution functions JTIJ, 

fs(q s ,Ps,r) = — ^ / rfr s F.^g^e*-, (2.11) 

F s (q' s ,qs,r) = J dp s f s (q s , Ps , T)e~^, (2.12) 

where 

1 „ 1 

Qs = qs- -r„ q s =q s + -r„ 

and g s = (<?i,g 2 , • • ■,?«)> = (pi,P2, • • • ,p s ), P s -r a = Y, s i= iPiPi^ dr s = dr x dr 2 ■ --dr s . 
One may substitute (2.12) into (2.10) and obtain the quantum BBGKY hierarchy of 
the Wigner distribution functions f s = f s (x s ,r), x s = (q s ,p s ), as 

± z£ [dx s+1 {,'"-■■ "- f 2 )/ ) ,,. :./;.! = 0. (2.13) 

.7=1 J 

Here dx s = dq s dp s = d 4 q x ■ ■ ■d A q s d i pi ■ ■ -d 4 p S) and the operators 9^ and 9j >a +i are 
represented as follows, 

9i- = ^(—-—), s+ i = ^—. (2.14) 
13 dq { \dpi <%/' dqi dp/ 
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For s = 1 and 2, one finds 



dfi 
dr 



dr 



Pi dfi 

+ T77T" + I 



dx 2 (e 



h,2/2 _ -26*1,2/2 



M dqi 

± if dx 2 (e i6 ^ 2 - e~^' 2 ) P lj2 f 2 = 0, 



h 



(2.15) 



+ E 



Pj dfi 
mi M % 



2 



3=1 
2 



0. 



(2.16) 



Equations (2.15) and (2.16) are exact. Since f 2 depends on f 3 , accurate solution of the 
hierarchy is impossible. One has, therefore, to apply some approximated approach. 
One of such approaches is the Bogoliubov one |]TB| , which we shall apply in the present 
consideration. 

According to the Bogoliubov hypotheses [[ITJ , 

1) It is possibe to find a kinetic state of any non-equilibrium system, provided 
that the average interval between two subsequent collisions is much longer than the 
duration of the collision. In this kinetic state, 



f s (xt, ...,x s ;t) = f s (x!, ...,x s \fi), 
dfi 



dr 



(2.17) 
(2.18) 



2) There are no correlations in the initial state of a system. One can introduce 
the displacement operator, 



V°f s {x° l ,...,x° s )=f s {x 1 ,...,x s ), 



(2.19) 



where the values of each x at r = 0, and their values at 

r. The non-correlative condition at the initial state implies 



' —T 



f s (x u ...,x s )~ n fi(xj) 

l<j<s 



0. 



(2.20) 



Starting from the Bogoliubov hypotheses, it is possible to derive a kinetic equation. 

Although the invariant interaction potential has infinite support in space-time, 
since it depends on (q 1 — q 2 ) 2 — c 2 {ti — t 2 ) 2 , its long-range part is necessary close to 
the light cone. It has been shown p3| , that wave operators exist in scattering theory 
if the support of the wave function does not extend to zero mass. The space-time 
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volume v of the effective interaction is therefore bounded. We shall assume here 
that it may be taken to be small, as in the first hypothesis of Bogoliubov. One can, 
therefore, write 



f£=A°{x\f 1 )+vA 1 (x\f 1 ) 



(2.21) 
(2.22) 



L /? • <•./; • <-./; • • • • • 
In the first-order approximation, one sets 

/ 2 = / 2 ° = /i(l)/i(2) (2.23) 

(henceforth we use the notation 1 = (x\\ r), 2 = (x 2 \ t), etc.) and finds from (2.15) 



0/1(1) + pi0/i(i 



dx 2 (e* 1 - 



/2 _ e -ifl 1>2 /2 



/i(l)/i(2) 



<9r ' M d qi 

± i J dx 2 (e i6 ^ 2 - e^ 2 ) P 1 , 2 f 1 (l)f 1 (2) = 0. 



(2.24) 



This self-consistent equation is a relativistic generalization of the quantum Vlasov 



equation 15 



In the second-order approximation, one writes a formal solution, 
f,(x 1 ,...,x a \f 1 )= 9{x h Xj) f x {k), 

i<j<s k^i^j 

where 

is a two-body correlation function, whose boundary condition is 

lim V_lg{xi,Xj) = 0. 



(2.25) 



(2.26) 



(2.27) 



Eq. (2.25) means that s-body effects are correlated by two-body effects. One can 
write 



dh = dhdfi 
dr dfi dr 



A°{x\h)+vA l {x\h) 



D f$ + v\D g(x 1 ,x 2 ) + D l f$ 



(2.28) 



where 



dfi df x 
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One now uses Eqs. (2. 16), (2. 25) and obtains 
D g( Xl ,x 2 ) + j^^-JLgix^x*) + iY,{e**l*-e-**l*)g(x x ,x 2 ) 

3=1 M ° q i 3=1 

= - i(e<*' 2 - e-^/j /i(l)/i(2) - iJdx 3 (e<*' 2 -e-<*' 2 
xg(x 2 ,x 3 )f 1 (l) - ijdx 3 (e i9 y 2 -e-<^f 1 (2)g(x 1: x 3 ) 
T ijdx 3 (e<*' 2 - e -<»/ 2 ) / 1 (1)/ 1 (2)/ 1 (3) 

=F </dx 3 (e<3/2 _ e -<»/ 2 ) / 1 (1)/ 1 (2)/ 1 (3). (2.29) 
Once x 2 ) is known, one can obtain the two-order-approximated equation for f\ : 

+ + i(e^ /2 -e-^/ 2 )h{l) + i[dx 2 ( e <^ 2 -e-<^ 2 ) 

or M oq 1 v ' J \ J 

*9(xi,x 2 ) ± ijdx 2 [e<^ 2 -e~<^f l {l)f l {2) = 0. (2.30) 



Here 



and 



a' l a a' l a dU x d 

^1,2 - _y l,2, ^1,3 - - y l,3, ^1 - 



f v oqi opi 

1 



C/i(?i,r) = ~Jdx 2 h(2)V( qi -q 2 ) (2.31) 



is the mean-field potential. In general, it is very difficult to obtain simultaneously 
solutions of Eqs. (2.29) and (2.30). In the following section we show how Eq. (2.29) 
can be solved for a quasihomogeneous system. 

2.1 Quasihomogeneous system 

The condition on a quasihomogeneous system is 

g(xi,x 2 ) = g(qi - q2,Pi,P2) = g(q,Pi,P2), (2.32) 

i.e., the correlation function depends only on the relative coordinates. In this case, 
one obtains a formal solution for g(q,pi,p 2 ) by means of the displacement techniques, 
as follows: 



CO 



i ^ e 2a g ^a P1 dv2> — e 28 ? l *i )V\q — — — — r 



9(q,Pi,P2) = / dr 
Jo 

x/i(l)/i(2) + i J dq'dp 3 (e^'&i - e - ^ 3 *^ V (q - q 
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Pi -P2 



rj x (g(q,p 2 , P3 )h(l) ± A(1)A(2)A(3)) 

dqdpz \ e 29qdp 2 — e 2 s< ? a n \ V \q — q — — 

x (g(q',Pi,Ps)fi(2) ± / 1 (1)/ 1 (2)/ 1 (3)); 
= ^Jdr [( e < 2 / 2 -e-< 2 /^/ 1 (l)/ 1 (2) 

+ J dx 3 (e<a/2 _ e -<»/ 2 ) x (^(a; 2 ,xs)/i(l) ± / 1 (1)/ 1 (2)/ 1 (3)) 

± / rfx 3 (e<a/2 _ e -< 8 /*) x (^^3)^(2) ± / 1 (1)/ 1 (2)/ 1 (3)) 

(2.33) 

In order to solve Eqs. (2.30) and (2.33), it is convenient to introduce the Fourier 
transform, as follows: 

~g(k, Pl ,p 2 ) = Jdq g(q,p 1 ,p 2 )e- ik - q , (2.34) 

V(k) = Jdq V(q)e- ik - q . (2.35) 
Substituting Eqs. (2.34), (2.35) into (2.30), one finds 

f +Sf («**-«-•*) <w 

where 

= J dp 2 g(k,p 1 ,p 2 ), (2.37) 

F Wi = 1 ^ m/2 ~ e ~ im/2 ^ ± 1 1 dX2 ( e<2 " 6 ^ ,2 ) /i(1)/i(2) - (2 - 38) 

Making a Fourier transform of Eq. (2.33), one obtains, after some manipulations, 



im/#,Pi) = / d P2 J Vl 'f\ l2 [tfmrw-rtmrw 

J k\lTV 2 , 3 ip\ 2 1 

Pi -p 2 



x5(k.^). (2.39) 



Here 



(2.40) 



2 / 

(the second sign ± in (2.40) distinguishes between bosons and fermions), 

f(p±^j = e ± **/( p ), (2.41) 
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and 



dp 3 



k ■ 



P1—P2 ' 

. M - 



— ie 



A + (3)-/r(3) 



(2.42) 



Substituting (2.39) into (2.36), one finally obtains 



Or 



+ 



Pi dfi(l) 



M d qi 

dp 2 S [k 



+ F 



dfi(l) 



71 



X 



dpi 

Pi ~P2 



(27T) 4 



dk 



k d k d 

g 2 dp~i g 2 dpi 



/i + (l)/r(2) 



A + (2)/ 2 -(i: 



(2.43) 



Equation (2.43) is the kinetic equation of a gas of indistinguishable particles in the 
quasihomogeneous case (the improved BUU equation |15|). It reduces to the usual 
BUU equation provided that the many-body effects are neglected and that the first- 
order approximation for the term F is taken: 



dfi(l) 
Or 



+ 



Pi dfijl) 
M dqx 



dUidfijl) 
dqi dpi 



(2tt 



J dp 2 dp' 1 dp 2 5 4 (pi +p 2 - Pi - p 2 ) 



(piP2\Vi, 2 \p lP ' 2 )\ {/i(l')A(2') [1 ± /i(l)] [1 ± /i(2)] 
- / 1 (l)/ 1 (2)[l±/ 1 (l')] [l±/!(2')]}. 



^2.44) 



In contrast to the usual Boltzmann and BUU equations which are applicable in the 
restriction on the system to be dilute, Eq. (2.44) includes the influence of many-body 
effects. Therefore, Eq. (2.44) provides an essential improvement for the systems that 
have a higher particle density or a larger force range of particle interaction; e.g., for 
strongly interacting matter, heavy-ion collisions, or a cold relativistic plasma. 
Rewriting this equation in the form 



dm 

Or 



+ 



Pi dfi(l) 



dU x dh{l) 



M dqi 

x 

- /i(l)/i(2) 



(2*) 



^ J dp 2 dp[dp 2 5 A (pi +p 2 - pi - p 2 ) 



dqi dpi 

(pmlVi^p^ 2 {/i(l')/i(2) [1 + a fi(l)] [1 + ah{2)] 



l + o-/i(l') 



^/i(2')]} 



a 



±1, 



(2.45) 



one sees that it reduces to the usual Boltzmann equation for non-identicai particles 
for cr = 0. Thus, the three cases, 

1, Bose — Einstein, 
a = < —1, Fermi — Dirac, 

0, Maxwell — Boltzmann, 

can be treated by means of a unique equation, (2.45), which can be, therefore, called 
the generalized Boltzmann equation. 
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3 Boltzmann //-theorem 



We now wish to establish the relativistic analogue to the Boltzmann if-theorem 
and to prove that the entropy of an ensemble of events, evolving without external 
disturbances, is nondecreasing as a function of r. 

The density of states in phase space associated with the distribution n has been 
found in 

!{n + g — l)!/n!(<7 — 1)!, Bose — Einstein 
g\/n\(g — n)\, Fermi — Dirac 

g n /n\ Maxwell — Boltzmann 

where g is a number of states in each elementary cell of energy-momentum space 
(degeneracy) and n is the average occupation number. Assuming no degeneracy 
(g — > 1) and using Stirling's approximation 

IniV! « AHniV, N » 1, 

we obtain for the density of entropy in phase space, s, 



— = lnAr(n) 
Kb 



-nlnn + (1 + n) In (1 + n), Bose — Einstein 
-nlnn — (1 — n) In (1 — n), Fermi — Dirac 
-nlnn, Maxwell — Boltzmann 

(a + n) In (1 + an) — nlnn, o = ±1, 0. 



(3.1) 



Therefore, in the case we are considering, the entropy of the ensemble is defined by 
the functional 



S(r) 

kn 



dqdp (cr+ /i(g,p;r))ln (1 + cfi(q,p;r)) - /i(g,p;r) ln/i(g,p; r) . (3.2) 



Then, taking the derivative of S(r)/kB, using Eq. (2.45) and integration by parts of 
the space-time derivatives, we obtain, after some manipulations, 



1 dS 
ks dr 



71 



dqdp 1 dp 2 dp' 1 dp' 2 5 4 (p 1 + p 2 - pi - p 2 ) {piP2\V 1 ^ 2 \p l p 2 ) 



4(2tt) 12 

xfi(Q,Pi;r)f l (q,p2;r)f l (q,p' 1 ;T)f 1 (q,p' 2 ;r) 



1 



x 



+ a 



Ji(q,pi;r) 
l 



a 



x 



In 

- ln< 



JM^pwt) 

Ji(q,Pi,r) 



a 



a 



v /i(?>p'i;t) 

f >_ 

\fi(q,P2;r / 

( \ 

\fi{q,P2,T) 



Ji(q,P2; T ) 



+ o 



(3.3) 
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In the derivation of (3.3) the principle of microscopic irreversibility (e.g., detailed 
balance) 

/ / 2 ii 2 / / 

{V\V2\V\?\V\V<il dp x dp 2 = (Pi#jIVl,2|PiP2) dp 1 dp 2 
and the hermiticity condition 



were used. Since 



(PlP2|K,2|p'iP2> = (PlP2lK,2|piP2) 



/ / 2 i , 

{V\Vi\Y\,2\P\P2) <5 4 G°i +P2 — Pi — Pa) — 0' anc ^ the remaining factor 



in the integrand is non-negative, we obtain 

dS(r) 



dr 



> 0, 



(3.4) 



the relativistic Boltzmann if-theorem. 

This result implies that the entropy S(t) is monotonically increasing as a function 
of r, and hence the evolution of the system, as described by the generalized relativistic 
Boltzmann equation, is irreversible in r, but not necessarily in t. In a smooth average 
sense, one can argue that the entropy must increase in t as well. The support of the 
distribution function in t is finite at each r; as r increases, this supprort moves up the 
t-axis, since the system as a whole moves with the free motion of the center of mass. 
The entropy, according to the if -theorem in r, must therefore also be nondecreasing, 
in this coarse-grained sense, in t. In the nonrelativistic limit [24] t — > r, S(t) takes on 
the usual nonrelativistic form, and the nonrelativistic if -theorem for evolution in t is 
recovered. 

In the special case in which the ensemble consists of positive energy (or negative 
energy) states alone, a precise if -theorem can be proved for the Lyapunov function 



p r 

d 3 qdrl d A p — (cr+fi(q,p; r)) In (l+cr/i(g,p; r))-fi(q,p; r) ]n/i(g,p; r) 



ks J "* Jp°>o r M 
by the application of the arguments contained in ref. M. 



3.1 Relativistic four-momentum distributions 

As we have seen in the preceding section, the entropy (3.2) of a system of events 
increases, according to the generalized relativistic Boltzmann equation, monotonically 
in r. It means that the momentum distribution function monotonically approaches 
its equilibrium value f{ (q,p)- The equilibrium limit is achieved when 



dS(r) 
dr 



0. 



(3.5) 
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Since the integrand in (3.3) is definite, (3.5) requires that, for the equilibrium distri- 
bution fr(q,p), 

(3.6) 

this condition implies the vanishing of the collision term in the generalized relativistic 
Boltzmann equation (2.45). 

Since pi,P2 and p[, p' 2 are the initial and final four- momenta for any scattering 
process, the general solution of (3.6) is of the form 

ln ( 7(oT7 7 + ° J = Xi(v,P) + X2(q,p) + ■■■, (3.7) 

VA {q,p) ) 

where the Xi exhaust all quantities for which 

Xi(q,Pi) + Xt(l,P2) (3.8) 

are conserved in collisions. The quantities conserved in the sense of (3.8) are the indi- 
vidual event four-momentum p^ and mass squared m 2 = — p 2 (the latter is asymptoti- 
cally conserved in the scattering process [13]), and a constant (the one-event "angular 
momentum" M^ v = q^p u — q v p^ also satisfies this requirement, but does not change 
the structure of the result). Hence, the most general form of /} is given by [§, [| 

ln ( 1 +a]= -A(p - p c f -B, A = A(q), B = B(q), (3.9) 

VA (q,p) ) 

where p^Pcfi is an arbitrary linear combination of the components so that 

1 Bose — Einstein, 



fl°\q,p) 1 



e -A(p- Pc ) 2 -B _ a 



exp{-A(p-p c ) 2 - J B}-l' 



1 

cxp{-A(p-p c ) 2 -B}+l- 



Fermi — Dirac, 



e A(p- Pc f+B^ Maxwell - Boltzmann. 

(3.10) 

The physical properties of the distributions (3.10) are studied in || for the case of non- 
identical particles, and in || for the case of identical particles. We shall normalize 
these distributions as (the physical meaning of such a normalization is manifested 
below): 

Jdqdpfi°\q,p) = V^\ (3.11) 

where is the total four-volume occupied by the ensemble in space-time. Let us 
introduce the system of the space-time densities, as follows: 

dpfi%,p) = n?(q), (3.12) 
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dp x dp 2 f 2 i °\q 1 ,q 2 ,Pi,p 2 ) = n 2 °\q 1 ,q 2 ), (3.13) 

/#A*rfW™.*) - 4%,^), *, (3,4, 
Then the one-body density, nf\q), is normalized, in view of (3.11), as 

dqnf\q) = V ii) . (3.15) 

In the case of no g-dependence of A and B, Eq. (3.13) yields = 1. 

4 Mean-field potential. RMS 

In the equilibrium case, Eq. (2.31) for the mean-field potential entering the general- 
ized Boltzmann equation (2.45), reduces to 

U l (q 1 ) = i j dq 2 dp 2 fi°\q 2 ,p 2 )V( qi - <fc). (4.1) 

Averaging (4.1) over the ensemble gives, through (3.11)— (3.13), 

1 ( 

U = 7^7(iy J M\dpx /i (gi,Pi)l/i(gi) 

= 2VWv I dqidpidq2dp2 fi°\ ( l^Pi)fi > \ ( l2,P2)V(q 1 - q 2 ) 

= ^^Jdq 1 dq 2 nf{q 1: q 2 )V{q 1 -q 2 ) ) (4.2) 

where we have used the relation f 2 « /i(l)/i(2). 

The total energy density of the ensemble is defined by 

P = Po + Pint, (4.3) 

where po is the energy density of a free gas (no-interaction case) calculated in refs. 
||. [4|, and pi n t is the contribution of the interaction potential which is equal to 

Pint = N U, 

N being the particle number density per unit comoving "proper" three- volume V^> , 
Nq = N/V^. We now assume that for the interacting gas V' 4 ' /N ~ v; it then follows 
form (4.2) that 

f 

p = po+ 2^(4)^^(3) J d ^ dq2 n 2 (Qu Q2)V(q 1 - q 2 ). (4.4) 
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For a quasihomogeneous system, nf\qi, q 2 ) = nf\qi — q 2 ), so that (4.4) takes on the 
form 

p = po + 2V w v m J d i n 2\q)v(q). (4-5) 

By introducing hyperbolic variables for spacelike q, 

q° = q sinh /3, q 1 = q cosh /3 sin cos 0, 

g 2 = q cosh /3 sin sin 0, g 3 = gcosh/^cos^, (4.6) 
< q < oo, -oo < (3 < oo, < 6 1 < 7r, < < 27T, 
one can rewrite the integral in Eq. (4.5) as 

4tt j g 3 dgcosh 2 (3d(3 n { 2 ] (q)V (q) . 

This integral does not have, however, a simply interpretable nonrelativistic limit, as 



we discuss below after Eq. (4.12). Let us instead turn to ref. |19|, where the two- 
body relativistic quantum-mechanical bound-state problem has been studied. It was 
found that, if the support of the wave function of the relative motion is restricted 
to an 0(2, l)-invariant subregion of the full spacelike region, one finds a lower mass 
eigenvalue of the ground state than in the case when the support is in the full spacelike 
region. The solutions, moreover, have a simply interpretable nonrelativistic limit. 
This subregion was called by the authors the "restricted Minkowski space" (RMS). It 
has a parametrization (in contrast to (4.6) corresponding to the full spacelike region) 

q° = q sin 9 sinh f3, q 1 = q sin 9 cosh (3 cos 0, 

q 2 = q sin 9 cosh (3 sin 0, g 3 = q cos 9. (4.7) 

Clearly, q\ + ql~ql = q 2 sin 2 9 > (and q\ + q\ + q\ — g 2 = q 2 > as well). This sub- 
measure space is 0(2, l)-invariant, but not 0(3, l)-invariant. The representations of 
0(3, 1) are induced from the irreducible representations of 0(2, 1) which are provided 
by the eigenf unctions of the two-body bound-state problem |20[| . The fact that this 
restricted subregion admits a lower mass of the ground state than the full spacelike 
region constitutes a spontaneous symmetry breaking of the 0(3, 1) invariance of the 
dynamical equations. 

The restriction of the relative coordinates to the RMS corresponds to a restricted 
range of correlations available to the two events propagating in a bound state, i.e., to 
the range of q\—q 2 available at each r. In computing the full spectrum of the two-body 
problem the authors assumed that the wave functions of the excited states also lie 
in the 0(2, l)-invariant subregion, i.e., these correlations are maintained for excited 
states as well. Indeed, it was found that the partial wave expansion for scattering 
theory is recovered in this submeasure space as well |21|]. Here we shall assume that 
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this result has more generality and can be applied in statistical mechanics: for any 
two events, their mutual correlations lie in the relative 0(2, l)-invariant subregion of 
the full spacelike region. It then follows that the two-body density n 2 (qi — 92) will 
have support lying in the RMS associated with the relative motion gf — q^. Therefore, 
the integral in Eq. (4.5) will be nonzero only in the RMS associated with q, according 
to the nonvanishing support of the two-body density ^^(q). In this way we obtain 
for the integral in Eq. (4.5) 



J q 3 dq sin 2 6d6 cosh pd/3d<t>r$\q)V(q). (4. 



We shall also assume that, at any instant of r, the extent of the ensemble in the 
g°-direction is bounded so that = ■ At, where At is the range of the 
time variable for the system as a whole. Therefore, in Eq. (4.8) — ^ < q° < and 
integration on (3 gives 

Arcsinh(At/2(jsin6>) /\f 

cosh (3d (5 



-Arcsinh(At/2gsin6») q SU1 9 

Eq. (4.8) then reduces to 

At J q 2 dq sin 8d6d<j) n£\q)V(q) = AnAt J dq q 2 nf ] (q)V(q). (4.9) 

Using now the relations = ■ At and N/V^ = N , the particle number 
density, one finally obtains from (4. 5), (4. 9) 

/V 2 r 

P = Po + i L J d 3 qnf\q)V(q), (4.10) 
where d 3 q stands for 47rq 2 dq. In the same way it is possible to obtain the expression 



for the pressure of the interacting gas [18 



p = J)0 -^/A 9 ^4%). (4.11) 

We see that the expressions for p and p are precisely of the same form as those of 
nonrelativistic statistical mechanics [25], but with q = ^q^q^ replacing r = ^/q 2 ", 
and V(q) replacing V(r). The situation is quite similar to the one occuring in the 



two-body bound-state problem ||19|| , where, upon separation of variables in the RMS, 
one is left with a radial equation for q = y/q^q^ which is of the same form as a 
nonrelativistic radial Schrodonger equation for r = v^q 2 ". Separation of variables in 
the RMS therefore has a clear correspondence to the nonrelativistic problem, as first 
remarked in In the nonrelativistic limit, the relative variables q° and p° vanish 
(all the particles are synchronized in this limit ||), and the formulas (4. 10), (4. 11) 
acquire their standard nonrelativistic expressions p5*|. 
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We remark that the integral in Eq. (4.5) in the full spacelike region can be made 
convergent in the same way, by imposing the bounds on the time variable, as follows: 
—At/2 < gsinh/? < At/2. In this case integration on f3 results in the expression 



and similar relation for p. Hence, apart from Tav, the average passage interval in r 
for the events which pass through a small representative four-volume of the system 
||, contained in the expressions for po,po and N || there will be another T&y 
entering the expressions for p int and Pi nt , upon replacement for At in the correspond- 
ing formulas, through the relation (which represents the averaging of the equation of 
motion for q , dq°/dr = p°/M, over the ensemble, (E) being the average energy) 



In the nonrelativistic (or in the sharp-mass) limit, Tav — > oo, which provides a 
stationarity of the system in space-time, but not a non-trivial evolution in r . While 
p , p and N are preserved in this singular limit, due to the relation || T&yAm = 2n, 
where Am is the width of the mass deviation from its on-shell value, Pi nt and pi nt 
turn out to converge with T&y. Therefore, Eq. (4.12) and similar formula for p do 
not have a well-defined nonrelativistic limit, in contrast to (4. 10), (4. 11), which admit 
its clear form. This fact should be a source of a spontaneous symmetry breaking of 
the 0(3, 1) invariance in the correlation function of a many-body problem. 

We remark that no problem with the convergence of the integral in Eq. (4.5) 
arises in 1+2 dimensions (for the extent in the g°-direction bounded). Indeed, the 
3D analog of (4.5) reads 




and we obtain, in place of (4.10) 




(4.12) 




(4.13) 



By introducing hyperbolic variables for spacelike g, 



q = gsinh f3, q 



i 



q cosh (3 sin 9, q 2 = q cosh (3 cos 9, 



(4.14) 



< q < oo 
we rewrite the latter integral as 



-oo < 13 < oo, < 9 < 2tt 
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Integration on (3 gives 

/•Arcsinh(At/2g) /\f 

/ cosh (3d(3 = — ; 

J-Arcsinh(At/2g) q 

therefore, one obtains, via At = V {3) /V i2 \ N/V^ = N , 

P = Po + ^J d 2 qnf\q)V(q), (4.15) 
the 3D analog of (4.10) (d 2 q stands for 2nqdq), and a similar relation for p. 



5 Concluding remarks 

We have generalized the nonrelativistic approach of Yang and Yao, based on the 
Wigner distribution functions and the Bogoliubov hypotheses, to the relativistic case. 
We have derived the generalized Boltzmann equation which, in the case of indistin- 
guishable particles, improves the standard BUU equation in three main aspects: 

1) The effect of Pauli blocking, / — > /(l — / ), is included in the collision term. 
This is important for the collision processes at intermediate and low temperature, 
e.g., in heavy- ion collisions. 

2) The modified mean-field interaction is introduced into the collision term. This 
has a great influence on far-nonequilibrium states. 

3) The equation takes into account binary collisions corrected for many-body 
effects, wherein the many-body shielding effect can be obtained spontaneously. 

We have introduced two-body interactions, by means of the direct action poten- 
tial V(q), where q is an invariant distance in the Minkowski space. The two-body 
correlations are taken to have the support in a relative 0(2, l)-invariant subregion of 
the full spacelike region, in order to provide a good nonrelativistic limit to the basic 
thermodynamic quantities. Since the expressions for the energy density and the pres- 
sure are identical in form to those of the nonrelativistic theory, some of the results for 
the nonrelativistic interacting gas should be applicable for an interacting off-shell gas 
as well. For example, the equation of state of the ideal gas of non-identical particles 
is V — NqT; therefore, it follows from (4. 10), (4. 11) that the equation of state of a 
relativistic interacting gas should have the same form (by methods analogous to those 
of the standard cluster expansion [26]) as that of a similarly interacting nonrelativistic 



one, i.e. 
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N n T 



5>cn (A 3 iv c 



i=i 



where A = y2n/MT is the thermal wavelength, and ct/(T) is the Ith virial coefficient 
(^ = 1). 

Applications of the generalized Boltzmann equation to realistic physical systems, 
e.g., heavy- ion collisions, are now under consideration. 
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